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Elastic  scattering  by  a  spherical  inclusion 

Valeri  A.  Korneev  and  Lane  R.  Johnson 


Center  fur  Computational  Seismology,  Lawrence  Berkeley  Laboratory, 
and  Seismograph »c  Station,  University  of  California, 

Berkeley,  California  94720 

1.  Introduction 

The  earth  is  inhomogeneous  on  a  wide  range  of  scales  and  a  variety  of  methods  have  been 
developed  in  seismology  for  analyzing  the  effects  of  these  inhomogeneities.  The  theory  of  wave 
scattering,  as  developed  in  the  fields  of  optics  and  acoustics,  has  been  adapted  to  the  case  of  elastic 
waves  and  has  been  quite  useful  in  studying  certain  types  of  these  inhomogeneities.  For  instance,  Aki 
(1973)  used  scattering  theory  to  study  the  phase  and  amplitude  fluctuations  of  waves  arriving  at  a 
seismic  array,  Haddon  and  Cleary  (1974)  interpreted  the  precursors  to  PKIKP  as  due  to  scattering  near 
the  mantle-core  boundary.  Aki  (1969)  attributed  the  coda  waves  from  local  earthquakes  to  scattering  in 
the  lithosphere,  and  Aki  (1980)  considered  the  role  of  scattering  in  the  attenuation  of  waves.  In  parallel 
with  these  applications  of  scattering,  the  necessary  exteasions  in  the  theory  of  elastic  wave  scattering 
were  also  developed.  Korneev  and  Johnson  (1993a  1993b)  discuss  the  background  for  both  the  exact 
and  approximate  theoretical  developments  in  (his  area.  Of  particular  interest  to  the  subject  of  this  paper 
are  the  excellent  studies  by  Wu  and  Aki  (1985a,  1985b). 

The  full  treatment  of  elastic  wave  scattering  is  not  a  simple  task,  and  most  seismological  studies 
have  employed  various  approximations  in  their  ase  of  scattering  theory.  These  include  the  assumption 
of  only  one  type  of  wave  (acoustic  approximation),  the  assumption  of  a  low  contrast  in  material  proper¬ 
ties  (Born  approximation),  and  the  assumption  of  low  frequencies  (Rayleigh  approximation).  While 
these  approximations  appear  to  be  reasonable  in  many  cases,  a  rigorous  justification  of  their  use  is 
difficult.  One  method  of  checking  the  validity  of  the  approximations  is  to  compare  them  with  exact 
analytical  solutions.  The  purpose  of  this  paper  is  to  develop  and  discuss  the  properties  of  one  such 
solution,  the  scattering  of  plane  P  waves  and  S  waves  by  a  spherical  inclusion. 

A  spherical  inclusion  is  the  most  convenient  choice  as  a  test  model  for  comparison  with  approxi¬ 
mate  solutions.  It  is  one  of  the  few  objects  for  which  the  scattering  problem  has  an  exact  and  computa¬ 
tionally  tractable  solution,  and  it  has  the  desirable  property  of  being  describable  by  a  minimum  number 
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of  parameters.  The  treatment  of  the  canonical  scattering  problem  for  the  sphere  has  a  long  history.  For 
light  scattering  it  was  formulated  by  Mie  (1908)  in  terms  of  a  series  of  spherical  harmonics,  and  a 
comprehensive  discussion  of  this  topic  can  be  found  in  Van  der  Hulst  (1957).  Elastic  scattering  by 
spherical  obstacles  has  also  been  the  subject  of  many  publications,  with  some  authors  using  potentials  in 
their  approach  to  the  problem  (Ying  and  Truell,  1956;  Tmell  et.  al.,  1969;  Yamakawa,  1962;  Nigul  et 
al.,  1974;  Morochnik,  1983a,  1983b)  and  others  using  displacements  (Petrashen,  1946,  1950a,  1950b, 
1953;  Korneev  and  Petrashen,  1987).  The  present  paper  follows  this  latter  approach  and  a  detailed 
treatment  of  the  analytical  and  numerical  aspects  of  the  scattering  problem  for  P  waves  incident  upon  a 
spherical  inclusion  can  be  found  in  Komeev  and  Johnson  (1993a),  with  a  discussion  of  various  approxi¬ 
mate  solutions  in  Komeev  and  Johnson  (1993b).  These  results  are  extended  in  the  present  paper  to  the 
case  of  an  incident  S  wave  so  thru  comptirisons  can  be  made  between  the  relative  scattering  of  P  waves 
and  S  waves  by  various  types  of  spherical  inclusions. 

2.  Statement  of  the  problem 

Consider  a  two-put  isotropic  medium  consisting  of  a  spherically  symmetric  inclusion  Vx  (  part 
v  =  1  )  with  radius  r  =  R  having  elastic  parameters  A.t  =  X](r) ,  p,  =  p,(r)  and  density  p,  =  p,(r) 
which  is  embedded  in  a  homogeneous  elastic  surrounding  medium  (part  v  =  2  )  having  elastic  parame¬ 
ters  \  =  X1 ,  p  =  p2  and  density  p  =  p2.  The  inclusion  Vx  may  contain  a  number  of  internal  shells 
which  are  bounded  by  spherical  interfaces  where  the  material  properties  or  their  spatial  derivatives  are 
radially  discontinuous.  The  boundary  conditions  on  such  interfaces  as  well  as  those  at  the  surface 
r  =  R  are  linear  and  homogeneous.  We  assume  (hat  all  elastic  displacement  fields  under  consideration 
have  harmonic  time  dependence  of  the  form  e'm  where  to  is  the  angular  frequency.  Joint  Cartesian 
(x,y,z )  and  spherical  ( r.0,0 J  coordinate  systems  with  the  origin  at  the  center  of  the  inclusion  will  be 
used. 

Incident  from  medium  v  =  2  is  a  harmonic  disturbance  with  a  displacement  field  given  by 

£l0  =  U0  (x,y,z)  e"M  (2.1) 

The  interaction  of  (his  incident  wave  with  the  inclusion  gives  rise  to  additional  displacement  fields  both 
inside  and  outside  the  inclusion,  and  these  are  demited  by 

Uv  =  Uv(x,y,z)  ,  (v  =  1,2)  (2.2) 

Since  we  will  be  primarily  interested  in  the  properties  of  the  additional  disturbance  outside  the  inclu¬ 
sion,  this  field  with  subscript  2  will  be  referred  to  as  the  scattered  field  ■  U2 .  Thus,  the  total  field 


2 


U  in  the  outer  medium  v  =  2  Ls  a  sum  of  the  incident  wave  and  scattered  field 


U  =  Uo  +  U*  (2.3) 

Hie  field  U,  as  well  as  both  of  its  individual  components,  mast  satisfy  tbe  equation  of  motion  far  a 
homogeneous  isotropic  elastic  medium. 

(X+2p)  V2U  -  p  VxVxU  +  po)2U  =  0  (2.4) 

The  equation  of  motion  in  any  spherical  shell  within  the  inclusion  has  tbe  form 

(X,+2p,)  VJU,  -  p,  VxVxU,  +  ^  V-  Uj  t 

dr 


+  2-^^-  +  x  Vx  U.l  +  p.to’U,  =  0 


(2.5) 


We  denote  the  velocities  of  the  compressional  and  shear  waves  by 


'V"  *  ^ 


(2.6) 


We  require  that  the  scattered  field  satisfy  a  radiation  condition  at  large  distances  from  the  inclusion 


V,c  = 


-it  r  A,(0,$)  , 

~  •  «  p  >  -  -  ■  d*  * 


(r 


(2.7) 


where  kp  =  to/V/),2)  and  k,  =  to/V*21  .  The  functions  Ap(0,0)  and  A,  (0,q>)  will  be  referred  to  as  scatter¬ 
ing  diagrams  of  compressional  and  shear  waves,  respectively. 


3.  Spherical  vectors 

The  solution  will  he  developed  using  the  spherical  vector  system  of  Petrashen  (1945,  1949).  A 
fairly  complete  description  of  this  system  can  be  found  in  Korneev  and  Johnson  (1993a),  so  only  the 
essential  elements  of  the  system  will  be  listed  here.  Tbe  basis  vectors  for  the  system  are 

Y t  -  YjSUe,*)  =  r  x  VK,m(0,<|>) 

Y£  -  Y£(0,0)  =  (/+D  ?  r,m  (0.0)  -  r  VKfc,  (0,0)  (3.1) 

Y*  -  Y^(0.0)  =  /  f  Y,m  (0,0)  +  r  VY,„Am 
with  the  usual  definition  of  the  spherical  harmonic  functions 

Yb, (0,0)  =  <-“♦  PncosB)  .  /SO,  (-/  SwSl) 

The  vectors  of  this  system  are  linearly  independent  at  any  point  (0,0)  on  a  spherical  surface.  For  1=0 
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only  the  one  vector  Yd,  ■  f  is  nonzero. 


In  the  space  of  vector  functions  f(8,+)  defined  on  the  spherical  surface  ft 
OS0Sn  0  £  +  £  2x  dQ  =  sinG  dB  d+ 


the  basis  vectors  satisfy  the  orthogonality  relation 

|  y*£>  v/X  da  =  [cdf’]'2 

where  the  normalizing  coefficients  are  given  by  the  expressions 

J  2/  +  1 

*"  V  4x  /(/  +  1)  (/  +  m)! 

*"  V  4n  (I  +  1)  (/  +  m)! 

c~- 

*■"  V  4x/  (/  +  /«)! 

For  vector  functions  f(0,$)  with  a  finite  norm 


(3  2) 


(3.3) 


|  Ifl2  dCl  m  |  f*  f  dSl  <  « 

the  system  of  spherical  vectors  (3.1)  is  complete  in  the  sense  of  convergence  in  the  mean  for  a  general¬ 
ized  Fourier  series  expansion  of  f(9,$) 

-  i 

«e,*)  =  X  X  X  a,'*'  Y »*)  (3.4) 

«5*0.+,-  /■()  ma-V 

where 


=  [t-tiT']3  |  Y*/*’  •  f  c/n  (3.5) 

Using  the  completeness  of  the  vector  system  (3.1),  we  can  seek  a  solution  of  our  scattering  prob¬ 
lem  in  a  form  of  a  series 

U(r.e,<»=  £  d£?  Vi*\r)  YffW*)  (3.6) 

*j  jii 

Because  of  the  spherical  symmetry  of  the  present  problem,  the  3-D  scattering  problem  is  reduced  to  a 
1-D  boundary  problem  which  must  be  solved  for  the  radial  functions  yfH?(r) .  If  the  field  U  is  known 
on  any  spherical  surface  r  =  constant,  then  the  expansion  coefficients  of  (3.6)  can  be  determined  using 
the  orthogonality  of  the  spherical  vectors 
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Y*£*0*)  ■  U(r,0,#)  dCl 


(3.7) 


4.  Bask  expressions 

In  the  case  of  elastic  wave  propagation  in  a  medium  with  spherical  symmetry  a  critical  element  is 
the  traction  vector  on  a  surface  r  =  constant ,  which  has  the  form 

MU)  =  X  V-  U  f  +  2*1^-  +  Hi*  *  U]  (4.1) 

dr 

If  the  field  U  is  taken  to  have  the  form  (3.6),  then  the  corresponding  traction  vector  has  the  form 


MU)=X  WYfiKM) 

Kj  jk 


where  the  expansion  coefficients  are  given  by 
7 2(r)  =  d,{:  n 


3y£.  < 

dr  r 


7'«)(r)  =  d**.  i 


ft  +  2(i) 


3yt«i  Voo 
3r  +  r 


-  4(1 


7i:(r)  = 


21  +  1 


[(/  +  l)X  +  (3/  +  2)m] ~~  +  (/  +  2)[(/  +  1)X  - 


21  +  1 


^  -  a- 1)**. 

dr  r 


(/  +  1)(X  +  (1) 

’ayjj, 

21  +  1 

.  dr 

+  U  +  2) 


vi. 


diZ, 


21  +  1 


(/X  +  (3/  +  +  (/  -  !)[(/  +  l)(t  -  /x]-^ 


(42) 


(4.3) 

(4.4a) 


(4.4b) 


(4.5) 


Note  that  the  coefficient  d£!1  remains  coupled  with  the  same  vector  y£!’(8,$)  in  the  expressions  for  both 
the  displacement  (3.6)  and  the  traction  (4.2).  Differential  equations  for  the  radial  functions  y^*(r)  can 
be  obtained  by  substituting  (3.6)  into  (2.5)  and  using  the  orthogonality  of  the  spherical  vectors.  Solu¬ 
tions  in  the  form  of  power  series  for  the  general  case  can  be  found  in  Korneev  (1983),  but  in  this  paper 
the  emphasis  will  be  on  the  special  case  of  a  homogeneous  sphere. 

In  the  dynamic  theory  of  elasticity  it  is  useful  to  consider  displacement  field  as  a  sum  of  potential 
(P)  aid  solcnoida)  (S)  fields 
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U  =  V,  +  Ui 


(4.6) 


which  satisfy  the  conditions 

VxU,  sO,  V  •  Us  •  0  (4.7) 

and  represent  compression^]  and  shear  waves,  respectively.  Since  V  •  y(r  )Y£(0,4)  m  0  and 
V  x  y(r)Yoo  *  0  the  fields  (4.6)  have  the  form 


I 

120.1m  t 


FUr)  Yi  +/£<r)  Y £, 


/t>)Y£+/£(r)Y£,+/*(r)  Y£. 


4 


/il.UilS/ 

where  for  I  £  1  the  radial  functions  must  satisfy  the  equations 


(4.8) 

(4.9) 


*F* 

dr 


F,~ 

(/  -  1)—  - 
r 


(/ 


dF, *  F£, 

-  U  +  2)—  =  0  , 
dr  r 


+  (/  +  1) 


dft 

—•  +  (/  +  2)— — - 
or  r 


0, 


(4.10) 

(4.11) 


In  the  case  of  a  homogeneous  isotropic  elastic  medium  the  displacement  field  U  mast  satisfy  the 
equation  of  motion  (2.4).  Substituting  the  expressions  (4.8)  and  (4.9)  into  (2.4)  and  using  the  ortho¬ 
gonality  of  the  spherical  vectors,  one  obtains  dif  ferential  equations  of  the  second  order  for  the  radial 
functions.  These  equations  have  general  solutions  of  the  form 

f£(r)  =  a,l*jul(k,,r)  +  afchM(ker) 

Ft(r)  -  afcj,-x{kpr)  +  aQ;h,_,%r) 

fit(r)  =  hXjux(k,r)  +  b£h,.l{k,r)  (4.12) 

/ta(r)  =  l>i„ji-iikar )  +  bfchi.i(k,r) 
fl(r)  =  c£j,%r)  +  c£'h,(k,r) 

The  solutions  here  have  been  constructed  as  a  linear  combination  of  two  independent  solutions,  the 
spherical  Bessel  functions  y„ (kr )  and  the  spherical  Hankel  functions  of  the  second  kind  hm(kr).  Reids 
which  are  regular  at  the  origin  will  contain  only  the  spherical  Bessel  functions,  whereas  secondary  scat¬ 
tered  fields  which  must  satisfy  the  radiation  conditions  of  the  form  (2.7)  when  r  — ►  <»  will  contain  only 
the  Hankel  functions.  The  differential  equations  (4.10)  and  (4.11)  in  this  case  reduce  to 

=  -<C  .  (/  +  1VC  =  «C  .  (v  *  1,2)  (4.13) 
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We  assume  that  the  incident  wave  U0  is  regular  at  the  origin  so  the  radial  functions  of  this  wave 
will  contain  only  the  spherical  Bessel  functions.  Thus  the  general  case  for  the  incident  field  U0  is  given 
by  the  expression 


Uo  =  2}  cS,j,(k,r) Yi  +  + 

/jw 


(4.14) 


with  arbitrary  coefficients  for  the  P  disturbance  and  coefficients  bS»  ,  c£  for  the  S  disturbance. 
Introducing  a  set  of  "canonical"  incident  waves  defined  as 

p*.  =J/+1(V)Y,:,<e.<|>)  -  >m<V> Y,-(04) 

SVlm  =  /;,♦,(*,  r)Y,;.(e.<|»)  +  (/+!)>,.,  (ksr)Yi,(B4)  (4.15) 

SH,„,  =y,ajr)Y^(0.4» 


we  can  represent  (4.14)  as  the  linear  combination 


Uo 


(4.16) 


Each  of  the  waves  of  (4.15)  satisfies  the  equation  of  motion  (2.4).  The  wave  P!m  Is  a  pure  compres- 
sional  wave  and  S V;„,  and  SH,„,  waves  are  both  shear  waves. 

Now  coasidcr  the  incidence  of  canonical  waves  of  the  form  (4.16)  on  the  inclusion  Vt  .  Substitu¬ 
tion  of  the  field  U|  having  the  form  (3.6)  into  the  equation  of  motion  (2.5)  leads  to  a  separate  set  of 
differential  equations  for  radial  functions  for  tiny  pair  of  indices  l  jn .  Moreover,  the  equation  for  deter¬ 
mining  y£,(r)  separates  from  those  for  (r)  and  (r).  Also  note  that  and  azimuth  index  m  is  not 
present  in  any  of  the  coefficients  of  (he  differential  equations.  The  boundary  conditions  on  the  surface 
r  =  R  of  the  inclusion  tire  required  to  be  linear  and  homogeneous.  For  a  welded  elastic-elastic  inter¬ 
face  they  have  the  form 

U,  =  U»  +  UM.  and  !<1,[ul]  =  tr(2)[u(,  +  U„]  (4.17) 

with  the  usual  modifications  for  elastic-fluid  and  elastic-free  interfaces.  Again,  because  of  the  ortho¬ 
gonality  of  vectors  (3.1),  separate  boundary  equations  may  be  obtained  for  any  pair  of  indices  /,m,  and 
in  the  present  problem  these  equatioas  do  not  depend  upon  (he  index  m .  The  canonical  field  Pj*  will 
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excite  in  the  medium  v  =  2  a  scattered  tickJ  of  the  form 

U£  =  +  +  (4.18) 

The  canonical  field  SV/(B  will  excite  in  the  medium  v  =  2  a  scattered  field 

Ufa  =  [a, s%*,(fc,r) +  »/“*, +1(*,r)]Yi  +  [-a^Ai-^r)  +  </+l)6/“*l.,(*,r)]vil  (4.19) 

and  the  field  SH^,  will  excite  in  the  medium  v  -  2  the  scattered  field 

VL"  =  c?h,(k,r)  Yfi,  (4.20) 

The  set  of  coefficients  af‘p  ,  afF  ,  bf' ,  ft/*  ,  cf  which  are  contained  in  these  expressitms  will  be  called 
the  canonical  scattering  coefficients  for  the  inclusion  V,.  They  may  be  found  as  solutions  of  linear 
systems  following  substitution  of  the  relevant  expressi  t  )  the  boundary  conditions  and  using  the 
orthogonality  of  the  spherical  vectors.  Analytical  formulas  for  these  canonical  scattering  coefficients  for 

the  the  case  of  a  homogeneous  elastic  spherical  inclusion,  as  well  as  for  the  special  cases  of  a  fluid 

filled  spherical  inclusion  and  a  spherical  cavity,  are  given  in  Appendix  A.  These  will  be  discussed  in 
more  detail  later. 

Once  the  canonical  scattering  coefficients  are  known,  an  incident  field  (4.16)  specified  by  the 
coefficients  a£,  ,  b£,  ,  c,",  will  generate  a  scattered  field  which  can  be  written  as 

u* :  =  Y.  \  c&j‘ht(k,r)Y}l 
*• 

+  [[«/,>/''’  +  +  l[a&hF  +  *£*/“]  W*.')]y£, 

+  £“{#/»■«/£/"  +  ft/m«/*f  ]ft/-i(*,,'')  +  (/+l)[^toi^/«V  +  b£,bj£ j ht r )j Y^  •  (4.21) 

d 

This  represents  the  complete  solution  for  the  scattered  field  from  a  spherical  inclusion  for  an  arbitrary 
incident  wave. 

The  field  U(  inside  of  the  sphere  will  have  the  general  form  of  (3.6)  and  will  be  linearly  depen¬ 
dent  upon  the  source  coefficients  ,  ft,",  ,  c,"(  .  For  the  special  case  of  a  homogeneous  isotropic 
material  inside  the  inclusion,  U,  hits  the  the  same  form  as  (4.21)  with  all  of  the  functions  hk  replaced 
by  the  corresponding  functions  jk  and  with  a  new  set  of  canonical  coefficients  for  the  inner  medium 
v  =  1  (r  <  R )  .  For  the  sake  of  completeness,  analytical  expressions  for  this  internal  set  of  canonical 
coefficients  for  the  cases  of  elastic  and  fluid  spheres  can  be  found  in  Appendix  B.  However,  throughout 
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(be  remainder  of  this  paper  only  the  scattered  held  outside  of  the  inclusion  will  be  considered. 

One  method  of  determining  the  coefficients  a£,,  b£ ,,  and  cC  of  the  incident  wave  is  to  use  (3.7) 
and  integrate  the  product  of  the  incident  field  U(,  with  the  corresponding  spherical  vector.  Thus,  few  the 
case  of  a  plane  P  wave  propitgating  in  the  direction  of  the  positive  z  axis 

U0  =  e~*'Z  i  (4.22) 


and 


-i— 0+1) 

«£  =  «  2  8. 


}>n  .0 


fC  =  o 


cC  =  0 


(4.23) 


For  an  incident  plane  S  waves  propagating  in  the  direction  of  the  positive  z  axis  and  polarized  along  the 
x  axis 

U(,  =  x  (4.24) 


and 

,  I  f  1  -»■£<*♦ 1) 

<  =  ()  .  1C  =  -8m,Je  2 

„  21+1  f.,,  lxK  s  ) 

(-im  ~  ~  2l(l+l)  ^  •— t  8m.ije 

For  a  point  pressure  source  located  at  the  point  R(,  =  (2«,0,0)  where  z()  >  R 

it  r7  P 


(4.25) 


(4.26) 


and 

<  =  h,  (f,,Z„)5„,  d  ,  1C  =  0  .  C  =  0  (4.27) 

The  scattered  field  (4.21),  expressed  in  terms  of  spherical  unit  vectors  (f,0,$),  for  the  cases  of  the 
incident  plane  P  wave  (4.22)  and  the  incident  phuie  S  wave  (4.24)  are  given  in  Appendix  C. 

The  convergence  of  the  series  (4.21)  depends  upon  the  observation  distance  r,  the  canonical 
scattering  coefficients,  and  the  coefficients  of  the  incident  field,  and  each  combination  of  these  variables 
may  require  a  special  investigation.  The  basic  problem  is  to  estimate  the  number  of  terms  that  should 
be  included  in  the  series  in  order  to  achieve  a  certain  level  of  accuracy.  One  general  guideline  is  that 
the  number  of  terms  which  me  necessary  in  order  to  represent  the  incident  wave  on  the  surface  of  the 
inclusion  at  the  desired  accuracy  is  a  good  estimate  of  the  number  of  terms  required  in  the  solution 
series.  Korneev  and  Johnson  (1993a)  considered  this  problem  for  the  scattering  of  a  plane  P  wave  and 
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showed  that  (he  necessary  number  of  terms  in  the  series  could  be  estimated  by  the  formula 

/„  =  — +  N  (428) 

where  R  is  the  radius  of  the  sphere  and  N  is  a  constant.  A  value  of  N  =  15  is  sufficient  to  give  an 
accuracy  of  1(T*. 

5.  Flow  of  the  scattered  energy 

A  useful  method  of  characterizing  the  scattering  by  an  object  is  to  calculate  the  energy  of  the 
scattered  waves  and  compare  it  to  the  energy  of  the  incident  wave.  Various  forms  of  this  ratio  between 
the  scattered  and  incident  energies  are  called  scattering  cross  sections.  The  energy  of  the  scattered 
waves  can  be  obtained  by  calculating  the  energy  llux  of  scattered  waves  through  a  surface  S  that  com¬ 
pletely  surrounds  the  object.  Noting  that  the  energy  llux  through  a  surface  element  ds  having  a  normal 
n  is  given  by  (U  -  t„[U])  and  that  the  energy  llux  averaged  over  one  period  is  co  Im(U  -  t„*[U]}/2,  then 
the  total  energy  llux  per  period  through  the  surface  S  is  given  by 

F  =  j  Imj  (U  •  C(U])  ds  (5.1) 

where  (*)  denotes  the  complex  conjugate. 

Substituting  the  total  held  (2.3)  into  (5.1)  and  assuming  conservation  of  energy  (no  energy 
absorption  by  the  material),  we  obtain 

F  =  F„  +2F,  =  Imj  (U„.  •  C(U„])  ds  +  a  Imj  (U,f  •  C[U„])  ds  =  0  (5.2) 

where  Fk  is  the  total  energy  flow  of  the  scattered  held  and  Fc  describes  the  energy  of  coherent  interac¬ 
tion  between  the  scattered  held  U*.  and  the  incident  held  U(,  .  Physically,  the  phenomena  of  scattering 
describes  the  conversion  of  part  of  the  energy  of  the  primary  incident  wave  into  the  energy  of  the 
secondary  scattered  waves.  Therefore,  alter  the  incident  wave  has  interacted  with  the  inclusion,  it 
should  have  lost  part  of  its  energy.  However,  the  form;d  solution  (2.3)  leaves  the  incident  wave  undis¬ 
turbed.  This  means  that  the  additional  held  U„.  of  (2.3)  must  include  both  the  change  in  the  primary 
wave  along  with  the  secondary  scattered  waves.  We  will  return  to  this  problem  later  when  considering 
the  scattering  cross-sections  of  clastic  spheres. 

To  calculate  the  energy  How  Fx  of  the  scattered  held  (4.21)  caused  by  the  incident  held  (4.16) 
we  need  expressions  for  the  tractions  associated  with  both  of  these  fields.  Since  total  energy  flow  does 
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not  depend  from  die  shape  of  the  surface  S ,  we  let  S  be  a  spherical  surface  of  radius  r ,  which  is  arbi¬ 
trary  so  long  as  the  inclusion  is  contained  inside  S.  The  traction  vector  MU*.)  of  the  held  !)„  on  this 
surface  has  the  form 


► 

,OU  -  Jjc&ifc,  YZ  +  [[a£ar  + 

+  utf  j  A,~  -  (/ + 1 ) +  b&b^  ]fir]vu  |  (5.3) 


C,  =  Jit,  /■//,_,  (it,  r )  -  (/  +2)  hi  (k,  r  )j 

Vl  m  it  r  f</+2)/„+1(V>Y 

B,+  \  Ul+2)  hM{k,r)Y 

B,~  J  =  r  [k,r  hl(k,r)  ~  2  1  (/-l)  /i/-i(it,r)  J 


The  traction  vector  for  the  incidence  held  (4.16)  can  be  obtained  from  (S.3)  by  setting  all  canonical 
scattering  coeflicienls  equal  to  one  and  by  substituting  for  all  spherical  Hankel  functions  the 
corresponding  Bessel  function  in  (5.4M5.6). 

After  making  all  of  the  necessary  substitutions  in  (5.2),  performing  the  integration  over  S,  and 
some  tedious  manipulations,  we  obtain 

F  =  Fp  +  Fs 

1  sc  1  xr  T  1  jr 


=  2n(X  +  2\i)V'^J2l  +  1)|4 2),  \a*'a'PP  + 


+  2na  +  2n)Vl,'2j2l  +  l)||  ^ ^  +  Of  f  (5.7a) 


which  can  also  be  expressed  us 


* 

Fk  =  -  2wtt  +  2\i)Vp^J2l  +  |fltS|2Re{a/,|,J  +  R e^Of*} 

t  jn 

+  /(/  +  Dy*  rRefcf]  +  l/^l2Re{6f}  +  Re(r«V5)  ' 

L<2/  +  D-  JJ 


(5.7b) 


This  is  an  exact  result.  The  part  F£.  corresponds  to  the  energy  How  of  the  scattered  P  waves  and  the 
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pan  Fic  corresponds  to  (he  energy  How  of  the  scattered  S  waves.  As  can  seen,  the  combined  incidence 
of  both  P  and  S  waves  on  the  inclusion  can  cause  constructive  or  destructive  interference  in  the  scat¬ 
tered  field. 

The  equation  (5.7)  must  be  true  for  any  set  of  coefficients  for  the  incident  wave  (4.16),  which 
leads  to  the  following  four  independent  relations  between  the  canonical  scattering  coefficients 


lc,5l2  = 

-Re(cf) 

(5.8) 

In, "I2  + 

/(/  +  Uy'l/;,wl2  =  -Re{o,w>) 

(5.9) 

\afp\2  + 

/(/  +  Ur’lhfl2  =  -/(/  +  l)r'R e{h“) 

(5.10) 

2 

+  /(/  +  urW")  =  -  m  +  Ur Vs 

(5.11) 

Multiplying  (5.1 1)  first  by  ajiP  and  then  by  and  eliminating  the  quantity  a?phf>s  leads  to  the  expres¬ 
sion 

rial'll  +  2n,"f  =  l2U  +  nYjV'f  |»  + 

Using  the  equivalences  (5.9)  :uid  (5.10),  this  equation  reduces  to 

=  r-a  +  i)Y|/vf  (5.i2) 

The  equivalences  (5.8M5.IO)  will  be  used  in  (he  next  section  in  formulating  optical  theorems.  The 
last  equivalence  (5. 12)  will  he  used  later  when  considering  the  relation  between  P  -+  S  and  S  -*  P 
scattering.  The  equivalences  (5.8M5.12)  are  also  useful  in  verifying  the  accuracy  of  numerical  calcula¬ 
tions. 

It  is  worth  noting  that  the  result  (5.7)  would  also  have  been  obtained  if  the  radial  functions  had 
been  reduced  to  their  f;ir  field  asymptotic  expressions  before  substituting  into  (5.2).  This  means  that  the 
net  energy  (lux  due  to  (he  nettr-field  terms  in  the  solution  is  zero.  However,  as  shown  by  Korneev  and 
Johnson  (1993a),  these  ncar-lield  terms  can  significantly  affect  the  displacement  field  formed  in  the 
vicinity  of  the  inclusion. 

6.  Scattering  cross-sections  and  optical  theorems 

Here  we  consider  the  special  cases  of  tut  incident  field  consisting  of  either  a  plane  P  wave  or  a 
plane  S  wave.  Earlier  we  obhtincd  the  coefficients  (4.23)  and  (4.25)  which  represent  these  waves  in 
terms  of  the  spherical  vectors  (3.1).  Now  we  introduce  the  scattering  cross-section  o  as  the  ratio 
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a 


(6.1) 


fjL 

Fu 

which  is  the  energy  Dow  F,c  of  the  scattered  field  normalized  by  the  energy  flow  F0  of  the  incident 
wave  per  unit  area  normal  to  the  direction  of  propagation. 

We  begin  with  the  case  of  the  incident  plane  P  wave  (4.23).  For  this  wave  the  energy  flow  per 
unit  area  of  the  incident  wave  is 


Fq  =  a  +  2*txvy 


Substituting  coefficients  (4.23)  into  (5.7)  we  get  the  scattering  cross-section 


Ft 


=  o"  +  o*v 


=  t?  S(2/  +  ,HKT+  w  + 


(6.2) 


(6.3a) 


=  -jj  £<2/  +  1)  Re{tt/’/')  (6.3b) 

kf  /»> 

On  the  other  hand,  putting  the  coefficients  (4.23)  into  the  expression  (4.21)  for  the  scattered  field  and 
using  the  asymptotic  representation 


h,(z) 


1  -i({  -  ^(/+l)) 
—e 


<z  >  /) 


for  the  spherical  llankel  functions,  we  obtain  for  6  =  0 

p  e~*r' 

U'(0)  =  Ap(0)-—^ — i 

where 

APm  =  Y  (2/  +  1)  a," 

kf  i» 

Comparing  (6.3)  suid  (6.5)  we  have  the  equation 
ap  =  lin{4/>((») 


(6.4) 


(6.5) 


(6.6) 


which  is  the  optical  theorem  for  an  incident  phuie  P  wave.  This  equation  establishes  a  connection 
between  the  scattering  cross-section  and  the  amplitude  of  (he  scattered  field  in  the  forward  direction. 
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For  an  incident  plane  S  wave  (4.25)  die  procedure  for  obtaining  an  optical  theorem  is  similar. 
The  energy  flow  per  unit  area  of  the  incident  S  wave  is 


n  = 

Then,  using  the  coefficients  (4.25)  in  (5.7)  we  have 


(6.7) 


SL  =  +  <r“ 


pS 

^0 


2w 

"  Si-i  2 h 


rV  ft  KM) 


=  ~~r  2>  +  l)  Re{6»vv  +  tf) 
k‘  u  t 

The  forward  scattered  shear  wave  in  the  tar  field  has  a  form 

-ik.r 


a. 


II I  ! 


(6.8b) 


U£(0)  =  4,(0)- 


(6.9) 


where 


o*» 

4,(0)  =  -i-  £<2/  +  I)  [bF  +  cj] 
1 21 


(6.10) 


Comparing  (6.8)  and  (6.10)  we  have 
0s  =  ~Im{4,(0)| 

*5 


(6.11) 


which  is  the  optical  theorem  for  an  incident  plane  S  wave.  Optical  theorems  such  as  this  and  (6.6)  can 
be  useful  in  studying  the  attenuation  of  waves  due  to  scattering. 


7.  Comparison  of  P  -»  S  and  S  -»  P  scattering 

The  equivalence  (5.12)  allows  one  to  compttre  scattering  of  converted  waves  for  the  same 
scatterer.  Applying  (5.12)  to  (6.3a)  and  (6.8a)  we  see  that  the  scattering  cross-sections  for  converted 
waves  are  connected  by  the  simple  relation 


This  equation  says  that  the  scattering  cross-section  of  P  -*  S  converted  waves  is  significantly  larger 
than  that  for  S  — »  P  convened  waves.  This  result  is  valid  for  any  spherically  symmetric  scatterer. 
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In  practice  one  deals  with  the  amplitudes  of  scattered  waves,  and  so  it  is  useful  to  also  estimate 
the  mean  intensity  (squared  amplitude)  of  the  held.  Using  the  far  held  approximation,  the  mean  inten¬ 
sity  of  the  scattered  held  for  an  incident  plane  P  wave  has  the  form 


I,  = 


4JI, 


=  7^7  2  <2 /  +  »‘  |«/"f  +  10  +  DY2  |*/Wj2} 


p"  +  o” 
4nr3  4ayr2 


Ipp  +  IpS 


02) 


where  the  integration  is  taken  over  a  spherical  surface  at  the  radius  r  >  R  in  the  far  held  zone.  Simi¬ 
larly,  for  an  incident  plane  S  wave  we  have 


*  -  ifh 


da 


(21  +  1) 


*>,*>*» 

Tk-r- 


iu 


/(/+!) 


tii 


m 


ik  5>*»(kM«f] 

** '  in 


I.SP  +  (xs 


(7.3) 


Using  the  equivalence  (5.12),  we  have  for  the  ratio  of  the  mean  conversion  intensities 

hs  _  1  an  2_  _  ?  (X  +  2u)2 

lSP  "  r  <3SI‘  y*  M2 


Thus,  if  the  comparison  is  made  between  the  intensities  of  the  waves  rather  than  the  energy  flux  of 
(7.1),  the  asymmetry  in  the  average  conversion  between  P  and  S  waves  by  scattering  is  even  larger. 
Note  that  this  is  a  general  result  (ha  h  Ids  for  ;dl  frequencies.  For  the  "typical"  seismic  situation  where 
y  =  1/V3,  the  ratio  (7.4)  is  equal  to  18.  For  soft  media  where  Y  i-s  even  smaller  in  a  relative  sense,  this 
ratio  could  be  significantly  larger. 

For  the  case  of  homogeneous  spherical  inclusion  (elastic,  fluid,  or  cavity),  an  even  stronger  result 
can  be  obtained  which  involves  no  spatial  averaging.  Fust  note  that  the  converted  far  field  for  the 
incident  plane  P  wave  (4.22)  is  easily  obtained  from  (C.I)  and  has  the  form 


_  r 

n  i  * 


dPficoxQ) «  e 

— — - 9  -  Aps(b),8)^—^ — 0 


(7.5) 


15 


where  4w(co,0)  is  a  scattering  coefficient  fur  con verted  P  -» S  waves.  Similarly,  for  the  incident 
plane  S  wave  (4.24)  the  converted  far  field  can  be  derived  from  (C.3) 

=  -  i  cos*  V  f  l  +  a?  ^-P,\cos9)t  =  Asf((oM)1-^-t  (7-6) 

*'*  +  1  i  KpT  r 

with  the  scattering  coefficient  4s,  (to, 0,*)  .  Now  for  the  case  of  a  homogeneous  inclusion  it  follows 
from  the  solutions  listed  in  Appendix  A  licit 

at"  =  Yj/(/  +  1)*/"  (7.7) 

Using  (7.7)  in  (7.5)  and  (7.6)  we  have  the  relation 

4s,(d),6.*)  =  —  y22cos*-4w((i),0)  (7.8) 

Thus,  in  the  far  field  the  scattering  coefficients  of  the  converted  waves  have  the  same  functional  depen¬ 
dence  on  frequency  to  and  .uigle  6  .  For  the  case  of  an  incident  P  wave  the  problem  has  axial  sym¬ 
metry  so  Aps  has  no  dependence  upon  0,  but  this  is  not  true  of  Asf. 

Aki  (1992)  arrives  at  a  results  similar  to  (7.4)  and  (7.8)  using  a  more  general  approach  involving 
the  reciprocal  theorem,  lie  considered  only  the  case  where  the  polarization  of  the  P  and  S  waves  was 
in  the  same  plane  (0  =  0).  and  thus  obtained  tut  equation  similar  to  (7.8)  without  the  cos*  term.  Conse¬ 
quently,  because  the  average  value  of  cos3*  is  1/2,  his  equation  for  the  squared  amplitudes  does  not 
contain  the  factor  of  2  found  in  (7.4).  The  approach  of  Aki  (1992)  is  extended  in  Appendix  D  to  con¬ 
sider  polarized  S  waves  and  it  is  shown  that  general  reciprocal  relations  can  be  established  in  the  far 
field  for  an  arbitrary  localized  scatterer.  However,  it  appears  that  results  such  as  (7.8),  which  involve 
total  amplitudes  of  the  P  and  S  waves  incident  from  the  same  direction,  can  only  be  established  for 
scatterers  with  a  high  degree  of  symmetry. 

8.  Homogeneous  sphere 

The  results  that  have  been  presented  up  to  this  point  are  valid  for  any  inclusion  that  has  spherical 
symmetry.  To  proceed  further  requires  that  solutions  be  obtained  for  the  canonical  scattering 
coefficients,  and  in  order  to  do  this  it  is  necessary  to  specify  the  internal  structure  of  the  inclusion. 
Here  we  consider  the  special  case  where  the  material  properties  of  the  inclusion  are  independent  of  the 
radial  coordinate,  in  which  case  it  is  possible  to  obtain  analytical  solutions  for  the  canonical  scattering 
coefficients. 

The  scattering  of  plane  P  waves  by  a  homogeneous  sphere  was  treated  in  our  previous  papers 
(Korneev  and  Johnson,  1992a,  1993b)  where  detailed  interpretations  of  the  scattered  fields  were 
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presented.  Thai  set  of  solutions  has  been  expanded  to  include  the  scattering  of  S  waves  and  the  com¬ 
plete  set  of  canonical  scattering  coefficients  for  a  homogeneous  spherical  inclusion  is  given  in  Appendix 
A.  In  addition  to  the  elastic  inclusion,  formulas  are  also  given  fur  the  special  cases  of  a  Sukl  inclusion 
and  an  empty  cavity. 


Consider  the  solutions  given  in  Appendix  A  in  the  limit  of  low  frequency.  Then  the  most 
significant  scattering  coefficients  are  given  hy  the  asymptotic  expressions 
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The  coefficients  cf  and  cf  are  obviously  much  smaller  than  the  others  and  may  be  neglected  at  low 
frequencies.  The  case  of  a  fluid  inclusion  is  easily  obtained  from  (8.1H8.5)  by  putting  |t|  =0.  The 
coefficients  (8.1H8.5)  may  also  be  considered  for  the  case  of  intriasic  attenuation  iaside  the  inclusion 
by  assuming  that  the  elastic  parameters  X,  and  |t|  have  complex  values.  In  this  case  the  coefficients  for 
I  =  1  depend  only  upon  the  dcasity  contrast  of  the  inclusion,  whereas  the  other  coefficients  in  the  limit 
of  large  intrinsic  attenuation  go  to  the  values 
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These  fonnulas  represent  the  ease  of  a  small  absorbing  inclusion  which  captures  all  of  the  energy  which 
crosses  its  boundary. 

It  is  convenient  to  describe  the  energy  scattered  by  a  homogeneous  spherical  inclusion  with  a 
non-dimensional  normalized  scattering  cross  section.  Dividing  the  o  from  (6.1)  by  the  area  of  the 
geometrical  shadow,  we  have 
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For  an  incident  P  wave  in  the  low-frequency  limit  we  have 
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Similarly,  for  an  incident  S  wave  we  gel 
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For  the  low-contrast  cstse,  where 
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the  expressions  (8.12)  and  (8.13)  can  be  simplified  to 
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For  the  scalar  low-contrast  case  the  normalized  scattering  cross-section  may  be  described  by  the 
simple  formula  (Van  der  lluLst,  1057) 
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and  where  the  V,v)  are  chosen  as  either  Vp'  or  V'vv!.  according  to  the  nature  of  the  incident  wave.  This 
result  can  be  explained  by  the  interference  of  the  incident  and  refracted  waves  propagating  in  the  for¬ 
ward  direction,  where  the  panunctcr  a  is  just  the  phase  difference  between  these  two  waves  in  the  far 
field.  Motochnik  (1983a,  1983b)  derived  this  same  expression  for  the  low-contrast  elastic  case.  More 
recently  (Korneev  and  Johnson,  1993b)  comptired  this  result  with  the  exact  solution  for  an  incident 
plane  P  wave  and  found  grxnJ  agreement  for  contrasts  of  about  40%,  except  at  very  low  frequencies. 

It  is  clear  that  formula  (8.17)  is  asymptotic  to  the  value  2  in  the  high-frequency  limit.  This  is  the 
result  of  the  manner  in  which  the  problem  was  formulated,  whereby,  as  mentioned  in  section  S.  the 
secondary  diffracted  field  contains  both  the  scattered  waves  and  any  modifications  of  the  primary 
incident  wave.  For  the  perfectly  absorbing  sphere,  in  which  case  there  will  not  be  any  waves  that  are 
actually  scattered,  the  secondary  field  must  have  a  value  sufficient  to  cancel  the  incident  wave  in 
the  shadow  and  the  normalized  scattering  cross-section  will  have  an  asymptotic  value  of  1. 
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9.  Low-frequency  scattered  fields 


Expressions  (8.1)-(8.5)  for  the  canonical  scattering  coefficients  of  a  homogeneous  elastic  sphere 
may  be  used  to  obtain  low -frequency  asymptotics  for  the  scattered  field  (4.21).  Thus,  for  the  incident  P 
wave  (4.22)  we  have  the  far  field  asymptotic  solution 
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and  for  the  incident  S  wave  (4.24)  the  far  field  asymptotic  solution  has  the  form 
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where  the  following  iM)tation  has  been  used 
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The  case  of  a  low  contrast  between  the  material  properties  of  (he  inclusion  and  the  surrounding 
material  (Bom  approximation)  is  defined  by  the  conditions 
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and  then  the  expressions  (9.1M9.6)  become  the  s:une  as  those  obtained  by  Gubemalis  et  al.  (1977a 
1977b).  For  an  incident  P  wave  these  arc 
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and  for  an  incident  S  wave  they  are 
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10.  Numerical  results 

In  order  to  illustrate  some  of  (he  properties  of  the  solutions  derived  in  this  paper,  numerical  cal¬ 
culations  were  performed  for  a  few  sample  problems  involving  a  homogeneous  spherical  inclusion.  Far 
the  first  set  of  examples  the  material  properties  of  the  surrounding  medium  were  chosen  to  be  represen¬ 
tative  of  a  typical  continental  oust: 

V™  =  6.0  km  Is  ,  V®  =  3.5  km  Is  .  p:  =  2.7  g  Icin' . 

Five  different  models  were  used  for  the  inclusion,  with  the  properties  chosen  to  represent  a  variety  of 
different  types  of  obstacles  that  might  be  encountered  in  the  earth.  The  elastic  parameters  for  these  five 
models  are  as  follows: 


model  1 

- 

vy‘>  =  7.5  km  is  , 

V,("  =  4.4  km  Is  ,  . 

pi  *  3.1  g/cm 3  , 

model  2 

- 

V'J1’  =  4.5  km  Is  , 

V,0)  =  2.6  km  Is  , 

p!  =  2.3  g/cm3 , 

model  3 

- 

Vj"  =  3.4  km  Is  , 

V,(1)  =  0.0  km  Is  , 

p,  =  2.7  g/cm3  , 

model  4 

- 

V'<"  =  1.4  km  Is  . 

V,0)  =  0.0  km  Is  , 

p,  =  1.0  g/cm3 , 

model  5 

- 

V-;"  =  0.0  km  Is  , 

V,(U  =  0.0  km  Is  , 

pi  =  0.0  g/cm3  , 

For  each  of  these  models  the  scattering  problem  was  solved  for  an  incident  plane  P  wave  and  aLso  for 
an  incident  plane  S  wave.  The  results  of  the  calculations  are  presented  by  plotting  the  normalized 
scattering  cross-sections  and  ojj  as  a  function  of  the  parameter  4  =  £2  =  to/?  /Vj2). 

Models  I  mid  2  simulate  high-velocity  and  low-velocity  inclusions,  respectively,  with  the 
difference  in  materia!  properties  being  about  20%  in  each  case.  Figures  1  (model  I )  and  2  (model  2) 
present  the  normalized  cross  sections  for  these  two  types  of  inclusions.  The  general  pattern  of  the  total 
scattered  field  in  these  cross  sections  is  described  by  an  increase  as  co4  at  low  frequencies  which  merges 
into  long  large  oscillations  about  a  constant  value  of  2.0  at  higher  frequencies.  These  long  oscillations 
are  caused  by  the  interference  between  the  waves  that  propagate  through  the  inclusion  and  those  that 
propagate  around  it  (Van  der  Hulst,  1957),  and  the  asymptotic  value  of  2.0,  as  discussed  earlier,  results 
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from  the  fact  that  the  scattered  field  detains  botli  the  waves  scattered  by  the  inclusion  and  the  distur¬ 
bance  of  the  primary  field.  Superimposed  on  this  pattern,  particularly  evident  in  Figure  2,  are  some 
short  small  amplitude  oscillations  caused  by  multiple  reflections  within  the  inclusion.  Because  the 
low-velocity  inclusion  tends  to  focus  energy  within  the  obstacle  much  more  than  the  high-velocity 
inclusion,  these  short  oscillations  are  much  more  pronounced  for  the  low-velocity  inclusion. 

Of  particular  interest  in  Figures  1  and  2  is  the  comparison  of  the  scattering  crms  sections  for 
incident  P  waves  and  S  waves  at  low  frequencies  where  the  wavelengths  are  larger  than  the  size  of  the 
scatterer.  In  this  frequency  range  the  P  — »  S  scattering  is  much  stronger  than  the  S  -»  P  scattering,  in 
agreement  with  the  results  derived  in  section  7.  For  an  incident  P  wave  the  energy  scattered  into  the  S 
field  can  exceed  (hat  scattered  into  the  P  field,  whereas  in  the  case  of  the  incident  S  wave  the  amount 
of  energy  scattered  into  the  P  tield  is  negligible  compared  to  that  scattered  into  the  S  field. 

Models  3  and  4  are  lluid  inclusions,  with  model  3  simulating  an  inclusion  of  molten  mck  and 
model  4  simulating  ;ui  inclusion  tilled  with  water.  The  scattering  cross  sections  for  these  fluid  inclu¬ 
sions  are  plotted  in  Figures  3  and  4  and  show  a  pattern  similar  U>  that  of  the  elastic  inclusions  except 
that  all  of  the  features  are  shifted  toward  lower  frequencies.  Because  of  this,  the  scattering  reaches 
significant  levels  at  rather  low  frequencies  where  the  size  of  the  inclusion  is  still  much  smaller  than  the 
wavelength  of  the  incident  wave.  The  observation  made  for  the  elastic  inclusions  that  the  P  -»  S 
scattering  is  much  stninger  than  the  .S'  ->  P  scattering  is  even  more  pronounced  for  the  fluid  inclusions, 
with  the  scattered  S  field  dominating  the  scattered  P  field  at  low  frequencies  regardless  of  whether  the 
incident  field  is  a  P  wave  or  S  wave.  In  addition,  for  the  case  of  the  incident  P  wave  the  scattered  S 
field  is  now  comparable  to  the  scattered  P  field  over  the  entire  frequency  range. 

For  the  case  of  the  water  tilled  inclusion  (Figure  4)  the  resonant  features  of  the  scattering  cross 
sections  are  particularly  conspicuous.  11k*  positions  of  the  resonance  peaks  correspond  to  the  real  parts 
of  the  complex  rixXs  of  the  determinant  (A.*))  contained  in  the  denominator  of  the  canonical  scattering 
coefficients.  Some  of  these  rixits  (including  the  first  one)  inay  be  obtained  by  letting  /  =  I  in  (A.9), 
which  leads  to  the  equation 
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The  first  few  roots  of  this  equation  are  £i  =  2.1  ,  5.9 , 9.2  ....  . 

In  model  5  the  inclusion  is  a  hollow  cavity  ;uid  the  scattering  cross  sections  are  shown  in  Figure 
5.  It  is  useful  to  think  of  this  mtxlcl  as  a  modification  of  the  water-filled  inclusion  of  model  4  in  which 
VJJ,)  and  pi  are  reduced  to  zero.  This  helps  explain  why  the  scattering  cross  sections  of  Figure  5  are 


22 


essentially  smoothed  versions  of  (hose  in  Figure  4,  with  the  main  differences  being  related  to  the  fact 
that  the  cavity  has  no  resonances  associated  with  the  scattered  field  within  the  inclusion.  The  fact  that 
the  general  patterns  of  the  scattering  cross  sections  in  Figures  4  and  5  are  similar  indicates  that  this  pat¬ 
tern  is  controlled  primarily  by  the  vanishing  of  the  shear  modulus  within  the  inclusion. 

It  is  worth  pointing  out  that  the  properties  of  the  material  surrounding  the  inclusion  and 
wavelengths  of  the  incident  waves  are  identical  in  the  Figures  1-5.  However,  is  is  clear  that  the  com¬ 
mon  features  of  the  scattering  cross  sections  are  found  at  rather  different  frequencies  for  the  different 
types  of  inclusions.  This  result  can  be  explained  if  one  describes  the  frequency  dependence  of  the 
scattering  cross  sections  in  terms  of  the  wavelength  of  the  scattered  field  rather  than  the  wavelength  of 
the  incident  wave.  Note  that  in  applying  this  reasoning,  the  wavelengths  of  the  scattered  fields  both 
inside  and  outside  the  inclusion  must  be  considered.  This  general  principle  explains  why  the  scattering 
cross  sections  of  the  S  field  is  always  shifted  toward  lower  frequencies  with  respect  to  those  of  the  P 
field  (compare  the  upper  and  lower  piuicls  in  Figures  1,  2,  and  3).  why  low-velocity  inclusions  have 
scattering  cross  sections  that  me  shifted  toward  low  frequencies  with  respect  to  those  of  high-velocity 
inclusions  (compare  Figures  1,  2,  3,  4,  and  5),  and  why  the  position  of  the  resonance  peaks  in  the 
scattering  cross  sections  depend  upon  the  velocity  within  the  inclusion  (compare  Figures  3  and  4). 

This  same  type  of  reasoning  about  the  wavelength  of  the  scattered  field  also  helps  explain  the 
general  result  that  the  P  — >  S  scattering  is  stronger  than  the  S  -»  P  scattering  at  low  frequencies. 
From  section  9  it  is  clem  that  the  low-frequency  scattering  energy  is  proportional  to  (R /wavelength)*. 
Such  a  result  favors  the  scattering  of  S  waves  because  of  their  shorter  wavelength.  Another  way  of  say¬ 
ing  this  is  that,  using  the  scale  of  wavelengths,  an  inclusion  appears  larger  to  an  S  wave  than  to  a  P 
wave  and  thus  it  is  scattered  more  intensively.  What  is  not  so  obvious  is  that  the  ratio  in  the  scattering 
intensities  for  the  converted  waves  should  be  independent  of  frequency  and  proportional  to  the  squared 
ratio  of  the  velocities.  I1  wever,  it  is  clem  in  Figures  1  -  5  that  the  shape  of  o^v  curve  is  always  ident¬ 
ical  to  the  corresponding  curve,  with  the  iunplitudcs  of  the  curves  scaled  according  to  (7.1). 

In  Section  8  it  was  pointed  out  that  the  case  of  an  inclusion  with  intrinsic  attenuation  can  be 
treated  by  assigning  complex  values  to  the  elastic  panuneters  within  the  inclusion.  Examples  of  the 
normalized  scattering  cross  sections  for  this  type-  of  an  inclusion  are  shown  in  Figures  6  and  7.  The 
attenuation  was  characterized  in  terms  of  the  quality  factor  ()j  where 
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The  calculations  were  performed  lor  the  low  frequency  case  where  k,,R  =  0.05  and  the  figures  show 
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how  the  cross  section  changes  as  the  attenuation  of  the  inclusion  is  increased.  Figure  6  shows  the 
results  for  the  high-velocity  inclusion  listed  earlier  as  model  1,  and  Figure  7  is  for  the  low-velocity 
inclusion  of  model  2.  The  results  are  quite  similar  for  the  two  types  of  inclusions.  For  small  attenua¬ 
tion  the  scattering  cross  sections  arc  just  the  low  frequency  values  shown  in  Figures  1  and  2.  As  the 
attenuation  is  increased  the  scattering  cross  sections  also  increase  and  approach  the  limiting  values  asso¬ 
ciated  with  the  coefficients  given  in  (8.9)  and  (8.10).  The  attenuation  affects  the  scattered  S  waves 
much  more  than  the  scattered  P  waves,  as  the  scattering  cross  sections  for  the  incident  S  wave  reaches 
values  which  are  about  3  limes  diosc  for  the  incident  P  wave.  Furthermore,  for  both  cases  of  an 
incident  P  wave  and  an  incident  S  wave  the  scattered  held  consists  almost  entirely  of  S  waves. 

The  basic  solutions  presented  in  this  paper  are  completely  general  in  that  they  can  be  applied 
over  the  entire  frequency  range  and  to  ;ui  inclusion  of  any  size.  For  instance,  (he  scattering  cross  sec¬ 
tions  of  Figure  3  can  be  used  to  provide  a  rough  estimate  of  scattering  by  the  earth's  fluid  core.  More 
appropriate  results  can  he  obtained  by  choosing  the  following  material  properties  to  represent  the  earth’s 
mantle  and  outer  core: 

Vf]  =11.3  km  Is  ,  V/-'  =  6.2  km  Is  .  p2  =  5.0  gm  Ian 3 

v;11  =  9.9  km  Is  ,  V/"  =  0.0  km  Is  ,  p,  =  6.0  gmlan3 

These  velocities  were  chosen  to  match  (lie  average  travel  limes  through  the  mantle  and  core,  and  the 
densities  were  chosen  to  match  the  contrast  in  acoustic  impedance  at  the  mantle-core  boundary.  The 
radius  of  the  core  was  taken  as  3482  km.  For  this  example  it  is  instructive  to  consider  the  complete 
solutions  to  the  scattering  problem  in  the  time  domain.  The  expansion  coefficients  for  the  incident  field 
of  (4.14)  were  chosen  to  represent  a  point  pressure  source  at  a  radius  of  6300  km,  and  then  the  scat¬ 
tered  field  of  (4.21)  was  evaluated  and  transformed  from  the  frequency  domain  to  the  time  domain. 
(Korneev  and  Johnson  (1993a)  show  how  the  solution  for  a  point  source  is  easily  obtained  from  the 
plane  wave  solutions.)  The  spectrum  of  the  pressure  at  the  source  was  flat  below  a  comer  frequency 
corresponding  to  a  period  of  30  sec.  The  toed  solutions,  including  both  the  incident  and  scattered 
fields,  are  shown  in  Figure  8  at  6-degree  angular  intervids  for  a  ntdius  of  6371  km.  Note  that  this  is  an 
example  of  high-frequency  scattering,  as  kr  R  has  a  value  of  730. 

There  are  a  variety  of  interesting  features  on  die  seismograms  of  Figure  8,  but  the  discussion  here 
will  concentrate  primarily  upon  some  of  the  diffraction  effects.  A  good  example  of  this  is  the  arrivals 
that  fill  in  the  gap  between  the  PcS  and  PKS  phases.  The  geometrical  ray  arrivals  for  the  PcS  wave 
end  at  a  distance  of  72  degrees  and  those  of  the  PKS  wave  begin  at  122  degrees,  but  in  Figure  8  this 
gap  is  completely  filled  by  diffracted  waves.  AnoUier  exiunple  is  the  P2KS  phase  which  ends  with  a 
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caustic  at  a  distance  oi'  175  degrees,  but  strung  diffracted  waves  extend  out  to  180  degrees  and  back  to 
less  than  150  degrees  where  (hey  merge  with  the  PKS  phase.  The  situation  is  actually  more  compli¬ 
cated  than  this,  as  the  P3KS  geometrical  arrival  extends  out  to  1  IS  degrees  and  then  is  continued  by 
diffracted  waves  that  merge  with  the  P2KS  wave  nettr  180  degrees,  while  the  P4KS  geometrical  arrival 
extends  to  SS  degrees  and  is  coiuinucd  by  diffracted  waves  near  100  degrees.  Thus  the  waves  F4KS, 
P3KS,  P2KS,  PKS  and  PcS  along  with  their  diffractions  all  eome  together  to  form  a  complex  and  con¬ 
tinuous  group  of  waves  that  appear  tat  these  seismograms  between  a  distance  of  60  degrees  and  time  of 
2300  seconds,  extend  (Nil  to  180  degrees  and  1600  seconds,  and  then  continue  back  to  a  distance  of  0 
degrees  and  700  seconds.  Increasing  die  amplitude  of  the  seismograms  would  allow  this  same  type  of 
pattern  to  be  extended  to  include  the  P5KS  wave  and  other  higher  order  core  waves  of  this  family.  The 
same  type  of  phenomenon  also  occurs  for  the  PKP  family  of  waves,  although  these  waves  are  of 
slightly  lower  amplitude  duui  die  PKS  waves  suid  thus  not  as  easily  observed  in  Figure  8. 

The  distances  mentioned  above  for  the  regions  of  geometrical  arrivals  and  diffracted  waves  will 
be  slighdy  different  in  (lie  retd  earth  because  of  Hie  radial  variation  in  velocity  in  the  mantle  and  core. 
The  seismograms  will  also  be  considerably  more  complicated  because  of  the  additional  waves  caused 
by  the  inner  core,  the  surface  of  the  etuih.  and  S  waves  generated  at  the  source  in  the  case  of  earth¬ 
quakes.  However,  the  relative  amplitudes  of  the  the  different  waves,  the  distortions  in  the  waveforms, 
and  the  interaction  between  the  geometrical  and  diffracted  arrivals  shown  in  Figure  8  should  be  gen¬ 
erally  applicable  to  long  period  waves  in  the  earth. 

There  is  one  other  feature  present  in  Figure  8  which  is  worth  mentioning.  On  the  radial  com¬ 
ponent  at  a  dishuice  of  1 80  degrees  tuid  at  a  time  of  about  2500  sec  there  is  just  discernible  a  king 
period  wave  (period  of  about  6(X)  sec),  This  is  an  interlace  wave  of  the  Slonely  or  Schulte  type  which 
travels  on  the  mantle-core  boundary  with  a  velocity  of  about  4.4  km/sec. 

11.  Discussion  and  conclusions 

The  primary  purpose  of  (his  paper  is  to  present  in  a  convenient  form  the  exact  solutions  for  the 
scattering  of  P  waves  and  S  waves  by  a  spherical  inclusion  and  to  point  (Hit  some  of  the  important  pro¬ 
perties  of  this  solution.  However,  it  is  also  worth  considering  whether  these  results  can  be  used  to 
make  some  general  inferences  about  llie  scattering  of  elastic  waves  in  the  earth.  In  doing  (his  the  first 
point  which  must  be  discussed  is  the  applicability  of  results  for  a  spherical  inclusion  to  the  situation  in 
the  earth  where  die  shape  of  the  inclusion  is  often  unknown,  but  most  likely  different  from  that  of  an 
exact  sphere.  I  lere  one  cun  appeal  to  die  tact  that  scattering  by  a  sphere  represents  a  canonical  problem 
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for  a  more  extended  class  of  objects  with  relatively  simple  and  smooth  boundaries,  and  thus  reason  that 
these  results  should  apply  in  tut  approximate  manner  to  a  wide  class  of  objects  having  these  properties. 
In  the  low-frequency  range  (Rayleigh  scattering)  an  even  stronger  argument  is  possible,  as  it  was  shown 
in  section  9  that  for  this  case  die  solutions  depend  only  upon  the  volume  of  the  inclusion  and  not  upon 
its  shape.  Iltus  the  low  frequency  results  presented  in  this  paper  should  he  applicable  to  small  3D 
inclusions  of  any  shape  in  the  earth. 

In  the  low-frequency  range  them  is  a  strong  asymmetry  in  the  relative  scattering  of  P  waves  and 
S  waves.  The  P  — »  .S'  scattering  is  generally  much  more  intensive  than  the  S  -*  P  scattering.  This  is 
explained  in  a  qualitative  sense  by  the  fact  that  (lie  inclusion  appears  to  be  larger  to  the  S  wave  because 
of  its  shorter  wavelength,  and  the  Diet  that  the  scattering  is  controlled  by  the  wavelength  of  the  scat¬ 
tered  wave  rather  lh;ui  die  wavelength  of  the  incident  wave.  It  is  common  for  an  incident  P  wave  to 
have  more  energy  in  the  scattered  S  field  thtui  the  scattered  P  field,  whereas  for  an  incident  S  wave 
almost  all  of  the  scattered  energy  is  in  the  S  lield.  This  suggests  that  the  coda  of  P  waves  should  con¬ 
tain  a  significant  proportion  of  S  waves,  while  the  coda  of  S  waves  should  be  predominantly  S  waves. 

This  asymmetry  in  the  scattering  conversion  of  P  stnd  S  waves  can  he  quantified  for  the  case  of 
the  spherically  symmetric  scatterers  considered  in  lliis  paper.  It  was  shown  in  section  7  that  the  mean 
intensity  of  the  P  —>  S  converted  waves  is  2  V*IV*  times  the  mean  intensity  of  the  S  -*  P  converted 
waves,  and  this  ratio  is  independent  of  frequency.  For  more  general  scafeters  it  is  possible  to  write 
reciprocal  relations  such  as  those  given  in  Appendix  D,  but  it  is  not  obvious  how  these  can  be  con¬ 
verted  to  intensity  ratios  such  as  that  just  given  for  a  spherical  scaltcrer.  However,  at  low  frequencies 
in  the  domain  of  Rayleigh  scattering  when;  only  the  volume  of  the  inclusion  is  important,  it  Is  conjec¬ 
tured  that  the  nitio  of  the  mean  intensities  of  the  converted  waves  will  approach  the  value  obtained  for 
spherical  scatterers.  Thus  this  strong  asymmetry  in  die  scattering  conversion  of  P  and  S  waves  is  likely 
to  be  a  general  result  when  the  wavelengths  arc  huge  compared  to  the  size  of  the  inclusion. 

The  scattering  from  it  lluid  inclusion  is  more  intensive  than  the  scattering  from  an  elastic  inclu¬ 
sion,  with  the  general  frequency  dependence  of  (lie  scattering  being  controlled  by  the  contrast  in  the 
shear  modulus.  Superimposed  upon  this  frequency  dependence  is  a  series  of  resonance  peaks  which  are 
controlled  by  compression^  velocity  of  the  lluid.  There  exists  the  potential  here  to  use  the  spectrum  of 
the  scattered  waves  to  estimate  the  dimensions  of  the  scatterers,  although  the  case  where  there  is  a  dis¬ 
tribution  in  the  size  of  the  scatterers  would  lend  to  smixxh  out  the  resonance  peaks.  Regardless,  the 
amount  of  energy  scattered  into  die  S  lield  by  an  incident  P  wave  is  an  effective  diagnostic  which  can 
be  used  over  the  entire  frequency  range  to  identify  lluid  inclusions. 
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In  the  ca.se  of  u  low  contrast  between  the  material  properties  of  the  inclusion  and  the  surrounding 
medium  only  a  few  of  the  scattering  coeflicients  need  he  included  in  the  low-frequency  range  and  they 
have  a  simple  dependence  upon  the  material  properties  (equations  (8.1H8.5)).  The  form  of  these 
coefficients  is  suitable  for  use  in  inverse  problems,  with  some  of  the  coefficients  depending  primarily 
upon  the  contrast  in  density,  others  depending  primarily  upon  the  contrast  in  shear  modulus,  and  others 
depending  upon  the  contrast  in  compressibility. 

One  feature  of  the  low -contrast  approximation,  as  is  true  of  most  Bom-type  approximations,  is 
that  it  does  not  satisfy  conservation  of  energy.  However,  this  Is  easily  remedied.  The  equations  (5.8)- 
(S.10)  are  derived  from  (5.2)  and  tire  essentially  statements  of  conservation  of  energy.  The  right-hand 
sides  of  these  equations,  which  involve  only  the  retd  parts  of  the  canonical  scattering  coefficients, 
represent  the  energy  icons  that  tire  coherent  with  the  primary  field  and  thus  account  for  the  change  in 
energy  of  the  prinuiry  field  dial  must  occur  when  additional  scattered  fields  are  generated.  Applying 
this  to  the  low-contrast  case,  we  note  that  the  coeflicients  (8.1)-(8.6)  are  completely  imaginary  and 
represent  only  the  scattered  fields,  die  retd  parts  which  represent  the  change  in  the  primary  field  having 
been  dropped  in  the  approximation.  However,  these  real  ptuts  can  he  recovered  from  the  expressions 
(5.8M5.10),  and  including  the  retd  parts  will  restore  the  conservation  of  energy.  Note  that  because  the 
coefficients  (8.1)-(8.6)  tdl  have  an  (oj)1  frequency  dependence,  the  real  parts  of  these  coefficients  will 
have  an  (to/’  frequency  dependence,  which  in  most  cases  will  make  them  small  enough  to  be  neglected. 
However,  in  some  situations,  such  as  studies  of  attenuation  of  primary  waves  due  to  scattering,  these 
real  parts  of  the  scattering  coefficients  should  be  included  in  order  to  achieve  a  formulation  more  com¬ 
patible  with  energy  conservation. 

A  caveat  involving  intrinsic  attenuation  should  be  mentioned  here.  As  mentioned  in  section  8, 
the  canonical  scattering  coefficients  ctut  be  modified  to  include  intrinsic  attenuation  by  introducing  com¬ 
plex  elastic  parameters.  However,  in  this  case  some  of  the  relations,  including  (5.8M5.11).  are  no 
longer  valid  because  strain  energy  is  no  longer  conserved.  The  analysis  of  the  low-contrast  approxima¬ 
tion  is  still  possible,  as  introducing  complex  elastic  parameters  into  (8.1M8.6)  produces  real  parts  of 
these  coefficients  which  are  proportioned  to  (to)'  and  which  cause  an  attenuation  of  the  primary  field 
due  to  the  inlriasic  attenuation  which  dominates  then  due  to  the  scattering. 
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Appendix  A 

Consider  the  case  of  a  homogeneous  elastic  sphere  of  radius  R  with  elastic  parameters  A,,  |it  and 
density  p(  surrounded  hy  a  medium  having  elastic  parameters  Xj,  p2.  and  density  p2  with  the  continu¬ 
ous  boundary  conditions  (4.17).  The  ctutonical  scattering  coefficients  have  the  following  analytical 
representations. 
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An  expression  for  A  pp  c;ui  he  derived  from  (A3)  by  substituting  for  the  functions 
M£2)  =  /-1././+I)  the  corresponding  functions  -yt(^2).  Analogously,  A-V.y  can  he  derived  from  the 

same  expression  by  substituting  for  the  functions  /i*(t|2)  (ft  =  /-U./+1)  the  corresponding  functions 
-A  0)2)-  The  following  definitions  have  been  used  in  equations  (A2)-(A5). 
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For  the  case  of  a  fluid  within  the  sphere  (p,  =  0),  the  above  expressions  reduce  to 
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Expressions  lor  and  for  A„  can  be  derived  from  A  in  the  same  way  as  in  the  elastic  case. 
For  the  case  of  a  cavity  the  above  expressions  can  be  further  simplified  to 
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Expressions  Air  A^>  aiul  for  A*v  can  he  derived  from  A  in  the  same  manner  ax  in  the  previous  cates. 


Appendix  B 

The  diffracted  field  inside  a  homogeneous  elastic  sphere  may  be  calculated  using  the  expression 
(4.21).  where  all  of  (he  spherical  llankcl  functions  ht  are  replaced  by  spherical  Bessel  functions  jk  and 


where  the  wavenumbers  kH  ;uhJ  k,  are  taken  for  the  inner  medium  v  =  1  .  The  canonical  scattering 
coefficients  for  the  inner  medium  have  the  loans 
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where  the  expression  for  A  is  given  by  (AT)  and  (he  notation  (A.6)  has  been  used. 
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For  the  case  of  a  lluid  within  the  sphere,  the  above  expressions  can  be  reduced  to 
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where  the  expressions  for  A  is  given  by  (A.‘J). 
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Appendix  C 

For  the  case  of  the  incident  plane  P  wave  (4.22)  the  scattered  held  in  the  outer  medium  has  the 

form 
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For  the  case  of  the  incident  plane  S  wave  (4.24)  the  scattered  held  in  the  outer  medium  has  the 

form 


in  terms  of  the  spherical  vectors  (3.1).  and 
V„  =  +  Vss 
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Appendix  D 

The  purpose  of  this  appendix  is  to  extend  the  analysis  of  Aki  (1982)  to  the  case  of  P  and  S  waves 
having  arbitrary  polarizations.  Consider  a  localized  scattcrer  with  a  spherical  coordinate  system  cen¬ 
tered  upon  it.  At  the  point  (r,8|,+|)  is  located  the  unit  radial  force 

t><r.e,.*,)  =  f  (D.l) 

The  scattered  S  wave  generated  by  this  force  and  observed  at  a  second  location  (r  ,83,43)  is 

•w (r ,63,^2)  —  Hni(r.e24i>  8  +  un2(r,B^2)  4  (D.2) 

=  uK(r  .83,43)  \(cos(a7)  8  +  sin  (o2)  4)) 

where  a2  is  the  polarization  angle  of  the  S  wave  at  the  second  location.  It  has  been  assumed  here  that 
the  distance  r  is  sufficiently  large  so  that  only  far  field  parts  of  the  solution  need  be  included.  At  this 
second  location  two  separate  forces  are  considered.  The  first  is  the  unit  transverse  force 

fstO-.02.42>  =  &  (D.3) 

At  the  first  location  (r  .8).+,)  this  gives  rise  to  a  scattered  P  wave  with  displacement 

•5ip(rl0,>+|)  =  uslP(r ,8, ,0| )  f  (D.4) 

The  reciprocal  theorem  suites  that 

fp(«’,8|,4i) '  —  f.v  iO  ,83,43)  ■  u^(r  ,63,43)  (D-5) 
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and  in  (be  present  case,  using  (D.1-D.4),  this  meaas  that 

Ujip(r  ,0|,ti)  =  u/%  i(r  ,02,^)  (D.6) 

The  second  force  to  be  considered  is  the  unit  transverse  force 

ts^r  ,02,^2)  =  ♦  (D.7) 

which  produces  a  P  wave  at  the  first  location  with  displacement 

■jv(r »8|,$|)  =  uxl/>(r ,6|,0|)  f  (D.8) 

Applying  the  reciprocal  theorem  in  this  case  yields 

usif(f  >®i>9i)  =  “w  2(r  .81  (D.9) 

The  two  expressions  (D.6)  and  (D.9)  are  the  recipnx:ul  relations  for  two  orthogonal  polarizations  of  the 
S  wave.  In  general  both  polarizations  will  he  present  in  a  scattered  P  wave  and  thus  both  reciprocal 
relations  are  required. 

In  the  vicinity  of  the  scatterer  the  amplitude  of  the  P  wave  incident  from  the  first  location  will  be 


ujf*  = 


4np  Vfr 


and  that  of  both  S  waves  incident  from  the  second  location  will  be 


4  4npV*  r 

Thus  (D.6)  gui  be  written  in  terms  of  amplitude  ratios  as 


uX\p(r  ,0|,$t)  Vx  a/..v,(r,02.<h) 


and  (D.9)  becomes 


uSlp(r  >01  >4*1 )  Kv*  .02,<^) 

Up  =  vj  p 


(Dll) 


These  are  exact  relationships  for  the  two  polarizations  of  the  S  wave  and  they  show  that  in  each  case 
the  P  -»  S  scattering  is  stronger  than  the  S  ->  P  scattering  by  a  factor  of  (VP/VS)2.  The  result  (D.10) 
with  =  0  is  essentially  that  derived  by  Aki  (1992).  While  these  results  are  very  general,  they 

are  not  entirely  suited  to  the  scattering  problem.  First,  they  deal  only  with  the  separate  components  of 
the  motion  and  cannot  be  converted  to  equivalent  expressions  involving  the  total  amplitude  of  the 
motion  unless  the  polarization  angle  of  uw  is  determined.  Guhcmatis  et  al.  (1977,  1979)  have  given 


general  formulas  for  the  far  field  scattered  P  and  S  waves  and  specific  formulas  for  a  few  special  cases 
of  homogeneous  inclusions  which  etui  he  used  to  determine  this  polarization  angle,  but  the  results  are 
not  particularly  simple.  The  second  problem  with  these  reciprocal  relations  is  that  they  involve  P 
waves  and  S  waves  incident  from  different  directions,  and  the  geometry  of  most  interest  in  scattering 
problems  involves  P  and  S  waves  incident  from  the  stune  direction. 

The  reciprocal  relationships  (D.10)  and  (D.ll)  can  be  further  simplified  in  the  case  of  symmetri¬ 
cal  scatterers.  In  the  case  where  the  scattering  object  possesses  cylindrical  symmetry  about  the  (0,,0i) 
direction,  the  polarization  angle  of  un  is  given  by 

sin  (0|  )sin  ) 


sm  (a-,)  = 


sin  (8) 


where 


cos  (8)  =  co.v(02)eos(0|)  +  .vi'fi(02)sin(0|>co.v(^2-^i) 

For  an  object  with  spherietd  .symmetry  this  result  holds  for  till  directions  and,  furthermore,  the  two  posi¬ 
tions  (0,.4>|)  and  (02,$|)  can  be  freely  intcrehtutged.  To  duplicate  the  problem  considered  in  this  paper, 
let  the  P  wave  be  incident  from  (0!  =  Ji,0|  =  0)  and  then 

UAV(r,6,0)  =  U/>.vl(r,0,0)  0  =  -  uHS(r,Q.$)  0 

Also  let  the  S  wave  be  incident  from  (his  stune  direction  (02  =  rt,4h  =  0)  and  without  loss  of  generality 
take  fj-2  =  0.  Then  use  (1X6)  tuid  the  spherical  symmetry  to  write 

uVf.(r.0,$)  =  Ups  |(r  ,k,0) 


=  mv(0)  upX(r,n.O) 


=  mv(0)  M.vi/.(r,0,^) 


Substituting  these  results  into  (1X10)  yields 


uv/.(r,0,0) 


u. s' 


<«> 


Kv3 

TTT  W 
ve 


Ups  ( r  ,0,0) 

Up' 


(D.  12) 


This  result  is  more  applicable  to  the  scattering  problem  as  it  involves  the  total  amplitude  of  the  P  and  S 
waves  and  both  the  incident  P  wave  ;uid  the  incident  S  wave  arrive  along  (he  same  direction.  It  is 
identical  to  (7.8)  which  was  obtained  from  the  exact  f:tr  field  solution  for  a  homogeneous  spherical 
scatterer.  Note  that  the  angle  0  in  this  result  is  just  the  tingle  between  the  polarization  vector  of  (he  P 
wave  and  the  plane  containing  the  polarization  vector  of  the  S  wave  and  the  scatterer. 
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P  wave  incident 


S  wave  incident 


Figure  1  Normalized  scattering  cross  sections  for  an  elastic  homogeneous  sphere  as  a  function  of  the 
parameter  kpR  =  (aR/Vp .  These  results  are  for  the  high-velocity  inclusion  which  is  listed  as  model  1  in 
the  text  The  top  two  panels  are  for  the  case  of  an  incident  P  wave,  while  the  bottom  two  panels  are 
for  an  incident  S  wave.  The  panels  on  the  right  are  expanded  versions  of  those  on  the  left  for  small 
values  of  the  argument.  The  dashed  line  represents  the  energy  scattered  as  P  waves,  the  dotted  tine 
represents  the  energy  scattered  as  S  waves,  and  the  solid  line  represents  the  total  scattered  energy. 
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Figure  2  Similar  to  Figure  1  for  the  low-velocity  inclusion  which  is  listed  as  model  2  in  the  text. 
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Figure  3  Similar  to  Figure  1  for  the  fluid  inclusion  which  is  listed  as  model  3  in  the  text 


Figure  6  Normalized  scattering  cross  sections  for  an  anelastic  homogeneous  sphere  as  a  function  of  the 
attenuation  quality  factor  {?*'  of  the  sphere.  These  results  are  for  the  high-velocity  inclusion  which  is 
listed  as  model  1  in  the  text  with  the  elastic  constants  within  the  sphere  modified  to  have  complex 
values.  The  frequency  Is  constant  with  kFR  =(oR/Vp  =  0.05.  The  top  two  panels  are  for  the  case  of 
an  incident  P  wave,  while  the  bottom  two  panels  are  for  an  incident  S  wave.  The  panels  on  the  right 
are  expanded  versions  of  those  on  the  left  for  small  values  of  the  argument.  The  dashed  line 
represents  the  energy  scattered  as  P  waves,  the  dotted  line  represents  the  energy  scattered  as  S  waves, 
and  the  solid  line  represents  the  total  scattered  energy. 
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P  wave  incident 


Figure  7  Similar  to  Figure  6  for  the  low-velocity  inclusion  which  is  Usted  as 
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Figure  8  Synthetic  seismograms  calculated  for  a  homogeneous  fluid  inclusion  that  simulates  the  earth's 
core.  The  seismograms  are  calculated  every  6  degrees  at  a  radius  of  6371  km.  The  source  is  a  point 
pressure  pulse  at  0  distance  and  a  radius  of  6300  km.  The  panels  on  the  left  are  the  radial  components 
of  motion  and  those  on  the  right  are  the  angular  components  of  motion.  The  upper  two  panels  are 
late-time  versions  of  those  below  with  the  amplitudes  increased  by  a  factor  of  10.  The  dotted  lines  are 
the  arrival  times  predicted  by  geometrical  ray  theory. 
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